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Abstract. A 3D simulation of the upper part of the solar convective zone is used to obtain information on the 
frequency component, Xfe, of the correlation product of the turbulent velocity field. This component plays an 
important role in the stochastic excitation of acoustic oscillations. A time analysis of the solar simulation shows 
that a gaussian function does not correctly reproduce the !/-dependency of Xfe inferred from the 3D simuation in 
the frequency range where the acoustic energy injected into the solar p modes is important ~ 2 — 4 mHz). 
The //-dependency of Xfc is fitted with different analytical functions which can then conveniently be used to 
compute the acoustic energy supply rate P injected into the solar radial oscillations. With constraints from a 3D 
simulation, adjustment of free parameters to solar data is no longer necessary and is not performed here. The 
result is compared with solar seismic data. Computed values of P obtained with the analytical function which 
fits best Xk are found ~ 2.7 times larger than those obtained with the gaussian model and reproduce better the 
solar seismic observations. This non-gaussian description also leads to a Reynolds stress contribution of the same 
order as the one arising from the advection of the turbulent fluctuations of entropy by the turbulent motions. 
Some discrepancy between observed and computed P values still exist at high frequency and possible causes for 
this discrepancy are discussed. 
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lent convective medium. The computed rates P at which 
the solar p modes are excited were found to be larger 
than those computed with a ID mixing-length solar model 
but still underestimate the solar seismic data by a factor 
~ 2.5. It was also found that the Reynolds tensor con- 
tributes about 20 % of the total acoustic energy injected 
into t he solar p modes, in cont rast with direct 3D estima- 
tions ijStein fc Nordlundll200lj) . These discrepancies were 
attributed to the assumed gaussian function for the dy- 
namic model of turbulence. 

In the present paper we therefore derive an empirical 
dynamic model of turbulence obtained from a 3D simula- 
tion of the upper part of the solar convective zone and then 
study the consequences of using this model on the com- 
puted excitation rates P. We compare our computation 
with solar seismic data and finally obtain on an improved 
model of stochastic excitation. 

The paper is organised as follows: The basic theoret- 
ical background and notations are recalled in Sect. 2. 
In Sect. 3, a 3D simulation of the upper part of the 
solar convective zone is used to characterise Xk in the 
domain where stochastic excitation takes place. The in- 
ferred i^-dependency of Xk is compared with the gaus- 
sian function and fitted with different non-gaussian func- 
tions. These functions are then used in Sect. 4 to compute 
the excitation rate P for radial p modes. The results are 
compared with solar seismic observations as provided by 
IChapl in et al.l l)l998j) and with computations in which the 
gaussian function is assumed. Sect. 5 is dedicated to dis- 
cussions and conclusions. 



2. Stochastic excitation 

2.1. The model of stochastic excitation 

We consider the model of stochastic excitation as de- 
scribe d in Paper I and assume here - as in Samad i et alJ 
^^2^)03^ - that injection of acoustic energy into the modes is 
isotropic and consider only radial p modes. Accordingly, 
the rate at which a given mode with frequency lvq is ex- 
cited can be written as : 
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In Eq. , po is the mean density, is the radial compo- 
nent of the fluid displacement adiabatic eigenfunction ^, / 
is the mode inertia (Eq. ifTOjl l. as = {dp/ds)^ where p de- 
notes the pressure and s the entropy, s is the rms value of 
the entropy fluctuations which are assumed to arise solely 
from turbulence, gr(^r)'') is a function involving the first 
and the second derivatives of w ith respect to r, <I> is a 
mean anisotropy factor defined bv lGoughl l)l977l) as 
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where u is the velocity field, < . > denotes horizontal 
average, () denotes time average, and w{r) is the mean 



vertical velocity {w = < u1 > — < u . 



Expressions 



for (7r(Cr: are given in I Samadi et al.l 11200: 

The driving sources ^/{(r, wq) and Ssir^oJo) arise from 
the Reynolds stress and the entropy fluctuations respec- 
tively : 

dfc E{k,r) E{k,r) 

k? Un Un 



SR{r,ujo) = 



Xki(^o,r) 



(3) 



dk E{k,r) Es{k,r) 
k^ Uq 

Xj 

dwXfc(wo + w, r)xfc(w,r) 



(4) 



where uo{r) = ^JWJiw is introduced for convenience, 
E{k, r) is the time averaged turbulent kinetic energy spec- 
trum, Es{k,r) is the time averaged turbulent spectrum 
associated with the entropy fluctuations and Xfc(a;,r) is 
the frequency-dependent part of the correlation product 
of the turbulent velocity field (see Section lT^ . In order to 
simplify the notation, we drop the explicit r dependence 
of the quantities in Eqs. (|1I4|I . 

2.2. The dynamic model of turbulence 

The dynamic model of turbulence is represented by Xfc('^)- 
In order to give a precise meaning to Xfe(w), we recall first 
some theoretical relations. Excitation by Reynolds stresses 
involves (fc, w), the Fourier transform of the second- 
order velocity correlations; here the indices i and j refer 
to any of the 3 directions of the velocity field. For incom- 
pressible, ho mogeneous and is otropic turbulence, 0ij(fc, lu) 
has the form lBatchelorlll97nl) : 



cj)ij{k,U!) 



E{k,u;) 
47rfc2 



(5) 



where E{k, lu) is the turbulent kinetic energy spectrum as 
a functio n of k and l o and 6ij is the Kronecker tensor. 
Following ISteiiJ lll96it . E{k,uj) is decomposed as 



E{k,uj) = E{k)xk{uj) 



(6) 



where Y k(io) satisfies the normalisation condition 
fTcnnckcs fc Lumlevlll98l Chap 8.1) : 



4-oc 



duJXkii^) = 1 



(7) 



According to the decomposition of Eq. (0 , Xk {^) is - at 
fixed fc - the frequency component of E{k,uj). 

According to Eq. ||SJ) and ©, Xfc(w) then represents 
the frequency dependence of (j)i^j{k,uj). In other words, 
Xk{'^) measures - in the frequency and fc wavenumber do- 
mains - the evolution of the velocity correlation between 
two distant points of the turbulent medium. 

The same decomposition of Eq. ^ is assumed for 
Eg{k,Lu). This leads to introducing xL the frequency- 
dependent part of the correlation product of the entropy 
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fluctuation. For simplifying the computation of P, as in 
Paper I, we assume xl = Xfc- We have checked that xl and 
Xk have almost the same behaviour in the region where 
excitation by entropy fluctuations is significant. 

In the present work, we consider only the excitation 
of radial p modes. Let Ez{k,uj) be the vertical compo- 
nent of the kinetic energy spectrum. We consider that 
Ez{k,oj) can be decomposed as E{k,uj) (Eq. Q). For 
isotropic turbulence we then have E{k,uj) — 3 Ez{k,uj), 
E(k) = 3Ez{k) and Xk — Xk, which is equivalent to stat- 
ing that the averaged and dynamic properties of the ve- 
locity field are the same in all 3 directions. 

The anisotropy factor $ - introduced in the expression 
for P, Eq. - partially takes into account the spatial 
and temporal anisotropy of the turbulence (<i> = 3 cor- 
responds to a isotro pic turbulence). It has been found in 
I Samadi et al.l l|2003f) that $ ~ 2 within the region where 
most of the excitation occurs. This shows that the time 
and space averaged properties of the medium are indeed 
anisotropic. One can therefore expect that the dynamic 
properties of turbulence differ between the horizontal and 
vertical directions. As the excitation of radial p modes 
is predominantly governed by turbulent elements moving 
in the vertical direction, an open question is whether one 
should consider xt rather than Xk in Eq- ® and Eq. 
when taking into account the dynamic properties of the 
turbulence. In the present work, we therefore characterise 
both Xk and Xfc from a 3D simulation and assess the con- 
sequences of using either xt or Xfc iri the calculation of 
P. 

2.3. A Gaussian function for Xk 

ISteinI l)l967|) and lMusielak et alJ l)l994l) suggested several 
analytical forms for Xki^j). The Gaussian Function (GF 
hereafter) is the simplest choice and is defined as 

Xfe(c.) = (8) 

where cok is its linewidth. 

In the time domain, the gaussian function, Eq. JH)), 
is the Fourier transform of a gaussian function whose 
linewidth is equal to 2Tk , where is a characteristic time 
correlation length. Hence Uk and Tk are related to each 
other as 
2 

tOk = — ■ (9) 

Tk 

The characteristic time Tk is usually associated with 
the characteristic correlation time-scale of an eddy with 
wavenumber k. As in .Balmforth (,1992^) . we define it as 

Tk^^ (10) 

kuk 

where the velocity Uk of the eddy with wave num ber k 
is rel ated to the kinetic energy spectrum E{k) by ijSteinl 
Il967^ 



r2k 



dkE{k) 



(11) 



The parameter A in Eq. (|10|l accounts for our lack of 
precise knowledge of the time correlation Tk under stellar 
conditions. 

I n the calculation of P. a GF is usual ly a s sumed for Xk 
fe.g. lGoldreich fc KeelevI Il977l) . BalmforthI lll992ft 'l. This 
assumption is equivalent to supposing that two distant 
points in the turbulent medium are uncorrelated. 

In Sect. 13 we use a 3D simulation of the upper part 
of the solar convective zone to derive the ^/-dependencies 
of Xfc and Xfc- Inferred i^-dependencies of Xfc and Xfe are 
compared to that of the GF. We next determine several 
analytical forms for Xfe and Xfe that can better represent 
their t/-dependencies. 

3. Constraints from the 3D simulation 

The analysis of a 3D simulation of the upper part of the 
solar convective zone provides constraints for several phys- 
ical parameters that enter the theoretical expression for 
the energy supply rate P injected into the solar p modes 
(Eq. The constraints may be considered to be of two 
types : 

• static constraints {static refers to spatial and time 
averages) determine the actual wavenumber depen- 
dency of E{k,z), the kinetic turbulent spectrum, and 
Es(k,z), the turbulent spectrum associated with the 
entropy. The static constraints also determine the 
depth profile of the wavenumber k^ at which con- 
vective energy is iniected into the t urbulent inertial 
range of E (as in lSamadi et alJ l|2003l) we assume that 
the wavenumber kg' , at which convective energy is in- 
jected into the turbulent inertial range of E^, is equal 
to fc^). They also provide the depth dependence of u^, 
the mean square velocity, w"^, the mean square vertical 
component of the velocity, s^, the mean square values 
of entropy fluctuations and $ = u^/w^, the rnean va l- 
ues of the anisotropy studied in lSamadi et alJ 1(20031) . 

• The dynamic constraints, on the other hand, concern 
the frequency component Xk and xl (^^^^ Sect. 12.2(1 . 

The s tatic constraints have been established in 
ISamadi et al. ( 200,'j) . Here we investigate the dynamic 
constraints. 



3.1. The 3D simulation 

We study a 3D simulation of the upper part of the so- 
lar convection zone obtained with the 3D numerical code 
developed by Stein & Nordlund (1998). 

The simulated domain is 3.2 Mm deep and its surface 
is 6 x 6 Mm^. The grid of mesh points is 256 x 256 x 163 
(i.e. ~ 23km x 23km x 37km), the total duration 27 mn 
and the sampling time 30s. 

Ou tputs of the simulation considered in lSamadi et alJ 
l|2003|) are the velocity field u(x,y,z, t) and the entropy 
s{x,y,z,t) where - as in lSamadi et al.l l '2003l) - z — t — Rq 
and Rq is the radius at the photosphere (i.e. where 
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T = Tefi). The quantities u(x, y, z, t) and s{x, y, z, t) were 
used to determine the quantities E{k, z), Es{k, z), w and 
involved in the theoretical expression for the excita- 
tion rate P. In the present work we use the velocity field 
y, z, t) to characterise Xk and Xk- 

3.2. Fourier transform 

At four different layers of the simulated domain, we com- 
pute the 3D Fourier transform, with respect to time and 
in the horizontal plane, of the velocity field u. These lay- 
ers cover a region where modes with frequency v ^2 mHz 
are predominantly excited. 

This provides z, v) where k is the wavenumber 
in the horizontal plane. Next we integrate ■u^(fe, z, v) over 
circles with radius k = \\k\\ at each given layer z. This 
yields u{k, z, v) and therefore E{k, v, z) = u^{k, v, z). The 
quantity Xk{v, z) is the frequency component of E{k, v, z), 
at fixed k. Hence, according to Eq.JHJ and 0, z) is 

obtained from v, z) as: 



Xk{v, z) 



E(k, V, z) 
J dv E{k, V, z) 



(12) 



where the integration over v is performed over the fre- 
quency range [— fmax, t'max] corresponding to the window 
of the Fourier analysis with respect to time (fmax — 
16 mHz). 

We proceed in the same manner for the vertical com- 
ponent of u. This then provides Ez{k,v,z) = u1{k,v,z) 
and Xk{v,z). 

3.3. Inferred properties of Xk snd xl 

Fig. n presents Xfc(^^) as it is obtained from the simulation 
for the wavenumber k at which E(k, z) peaks (fc = feg). We 
compare Xk{v) obtained from the 3D simulation analysis 
with the analytical functions GF, GEF and GLF (Eq. (jSJ, 
(O, and lO resp.). 

At the top of the superadiabatic region (for instance 
z = —0.4 Mm in Fig.l, this is the layer where the ex- 
citation is the largest), the GF does not correctly model 
Xk{v) (see Fig.nj. However the discrepancies between the 
GF and the simulation data occur mostly above the so- 
lar cut-off frequency {y ~ 5.5 mHz). Discrepancies be- 
tween the GF and the 3D simulation data have then mi- 
nor consequences for the p modes excitation in this region. 
This is not the case deeper in the simulation where the 
largest discrepancies between the GF and the simulation 
data occur in the frequency range where the dominant 
amount of acoustic energy is injected into the p modes 
(i/ - 2 - 4 niHz). 

To reproduce the shape of Xfc(^) obtained with the 
3D simulation, one needs a function which at high fre- 
quency decreases more slowly than the GF. For modeling 



the i^-dependency of Xk {v) , we thus propose three analyt- 
ical functions: the Lorentzian function (LF hereafter) 



Xk{i^) 



1 



7rwfe/2 1 + {2uj/ujk) 



(13) 



the Gaussian plus an Exponential function (GEF here- 
after) 



Xk{i^) 



1 



2 \ujk JH 



1 

2cji. 



(14) 



and the Gaussian plus a Lorentzian function (GLF here- 
after) 



Xk{^) 



1 



1 



.(15) 



All these functions satisfy the condition of normalisation 
of Eq. ©. 

We first assume a constant A = 1. As shown in Fig. ^ 
all these non-gaussian functions reproduce the z^- variation 
of Xk better than that obtained using a GF. 

In the middle of the excitation region (—0.5 Mm < 
z < 0.0 Mm) the overall best agreement is obtained with 
the LF. Below z ^ —0.5 Mm, the LF does not reproduce 
Xk well enough but still reproduces its z^-variation better 
than the other models. 

However we have so far assumed that A (or equivalently 
the eddy time correlation) is depth independent, which is 
a strong assumption. When A is allowed to vary with z, we 
find that decreasing the value of A below z < —0.5 Mm, 
the LF best models Xk below z ^ —0.5 Mm (e.g. A — 
1.6 at z=-0.64 Mm and A ~ 1.30 at z=-0.99 Mm, see 
Fig. . This shows that the variation with depth of the 
characteristic time (or equivalently the characteristic 
frequency tuk) is not correctly represented by the relations 
of Eqs. (|10lll(l when computed assuming a constant A 
below z ~ —0.5 Mm; increases faster with depth than 
expected from the relations Eqs. (|10lll(l . It is however 
found in Sect. 0] that this feature has negligible effect on 
P. 

The function Xk also decreases with the frequency 
more slowly than the GF (see Fig. Moreover, decreas- 
ing values of A for z < —0.5 Mm provide a better fit of 
Xk- But in contrast with Xk, Xk is overall better modeled 
with the GEF for z > -0.5 Mm and with the GLF for 
z < —0.5 Mm rather than with the LF (not shown). 

We conclude from the frequency analysis of the 3D 
simulation that the simple gaussian function cannot cor- 
rectly represent the actual dynamic properties of the tur- 
bulent medium. One may expect that the GF causes an 
underestimation of the acoustic energy injected into the 
solar p modes. Instead the frequency analysis favours a 
non-gaussian function for Xk and Xk that decreases more 
slowly with ly than the GF. 
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Fig. 1. The filled dots represent Xk{i^) obtained from the simulation for the wavenumber ko at which E{k,z) peaks. 
The solid curves represent the Lorentzian function (LF, Eq. the dots-dashed curves the Gaussian Lorentzian 
function (GLF, Eg. I15|l. the dashed curves the Gaussian Exponential function (GEF, Eg. I14|l and the long dashed 
curves the Gaussian function (GF, Eq.|HI). In these four analytical functions, A = 1 is assumed for the calculation oiujk 
fEgs. 1101 111) . Four different layers are considered: z — —0.04 Mm (the top of the superadiabatic region), z — —0.49 Mm, 
z = -0.99 Mm and z = -1.40 Mm. 



4. Consequences in terms of p modes excitation 

4.1. Computations of the excitation rate P 

Computation of the excitation rate P is performed as in 
ISamadi et al.1 1)2003|) except that here two analytical func- 
tions other than the GF are assumed for Xk, as discussed 
in Sect. 13 The computation process is summarised as fol- 
lows: The eigenfunctions (^r) and their frequencies (i/) are 
computed with Balmforth's (1992) non-adiabatic code for 
a solar ID mixing-length model based on Gough's (1977) 
non-local time-dependent formulation of convection. 

The quantities <&, w'^ and are obtained from the 3D 
simulation. The fc-dependency of E{k,z) is the Extended 
Kolmogorov Spectrum (EKS hereafter) defined as : 

E{k) OC (fc/fco) + ^ for fco > fc > /Cmin 

E{k) OC (k/ko)-^/^ for k>ka ^ ' 

In Eq. (|16|) . the wavenumber fco is the wavenumber at 
which E{k) peaks and fcmin is the minimal wavenumber 
reached by the 3D simulation (fcmin = 1-05 Mm^^). The 
variation with depth of fco is also given by the 3D simula- 
tion. The fc-dependency of the EKS reproduces the global 



features of E arising from the 3D simulation. The same 
model is considered for Es{k, z). E{k, z) and Es{k, z) sat- 
isfy the normalisation conditions: 



p + CC 

/ dfc£;(fc,z) = 1/2 $^ 

"'fcmin 

p + CC 

J dkEs{k,z) = 1/2 S2 



(17) 



The total energy contained in E{k,z) and Es{k,z), and 
their depth dependencies, are then obtained from the 3D 
simulation according to Eg. (|17|l . These theoretical esti- 
mates for P are then compared with the 'observed' P from 
IChanlin et all lll998l) 's seismic data, calculated according 
to the relation : 



I 



(18) 



where is the radius at which oscillations are measured. 



M 



dm ^j. 



(19) 
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Fig. 2. As in Fig. the filled dots represent z) obtained from the simulation at two different layers : z — 

—0.64 Mm (left panel) and z — —0.99 Mm (right panel). The other curves represent the LF fEa. 115(1 with different 
assumptions for A: The dashed curved correspond to A = 1 and the solid curves correspond to A = 0.80 a± z — —0.64 Mm 
(left panel) and A = 0.65 at 2 = —0.99 Mm (right panel). 



is the mode inertia and where the mode damping rate 
(7?) and the mode su rface velocity [vs] are obtained from 
IChaphn et aPlliggsl) . In Eq. lO, the mode mass I/^Kts) 
is given by the GMLT n iodel and we adopt r ,, = Rq + 
200 km consistently with lChaplin et all l|l998fl 's observa- 
tions. 

4.2. Comparisons with observations 

We investigate the effect of using different analytical func- 
tions for Xk (Sect. O in the computation of P . We first 
assume a constant value A = 1. Results are shown in Fig.O 
Computations performed with the GF underestimate the 
observed P values by a factor ^ 2.7. On the other hand, 
the LF, OFF and GLF choices result in larger values for 
the computed P than the GF one (~ 2 times larger) . This 
brings them closer to the observations, compared with the 
GF choice for Xk- The reason is that all the non-gaussian 
functions (the LF , the GEF and the GLF) - which in- 
deed better model x| a-nd Xk from the 3D simulation than 
does the GF - decrease more slowly with v than the GF in 
the frequency range where the mode amplitudes are large 
{y ~ 2—4 mHz). Consequently a larger amount of acoustic 
energy is injected into the modes with the non-gaussian 
functions than with the GF. 

In Sect. 13.31 the overall best models for Xk were ob- 
tained with the LF and with decreasing values of A below 
z ~ —0.5 Mm. We use a simple model for the depth vari- 
ation of A: 

A = 1 for z > -0.5 Mm 

A = 0.9 + 0.71(0.64Mm -I- z) for -0.5 > z > -1 Mm (20) 

A = 0.35 for z < -1 Mm 

We have computed P according to the simple model of 
Eq. I|20|l and assuming the LF. We find no significant 
changes for P compared to the calculations in which a 
constant value A = 1 is assumed (not shown). 



In Sect.Olwe found that Xk better modelled with the 
GEF for z > -0.5 Mm and with the GLF for z < -0.5 Mm 
rather than with the LF. However, as the stochastic exci- 
tation is the largest in the range —0.5 Mm < z < Mm, 
we can assume the GEF in all the domain. The LF results 
in a value for P,nax slightly larger than the one resulting 
from the GEF (only ~ 1.2 larger). A better agreement 
is then obtained with the analytical functions which fits 
best Xk (i.e. the LF) than the one which fits best xt (i-^- 
the GEF) in contrast with the intuitive idea mentioned 
in Sect. 12.21 that the excitation of radial p modes depends 
rather on the properties of Xk than on those of Xk ■ 

In the frequency range where observational constraints 
are available, differences between results obtained with the 
different adopted non-gaussian functions are of the same 
order as the actual error bars associated with the observed 
P values. 

Below the frequency range of the observations - i.e. 
below 1/ < 1.8 mHz - the differences between the dif- 
ferent non-gaussian functions are very large compared to 
the error bars (see bottom panel of Fig. OJ. Those dif- 
ferences are related to differences in the i/- variation of the 
non-gaussian models. Observational constraints at low fre- 
quency could therefore confirm that the LF is indeed the 
best representation for Xk{i^)- 

Important discrepancies still remain at high frequency 
{i' > 3.5 mHz). The excitation rate derived from the obser- 
vations decreases as ~ above v ~ 3.5 mHz) whereas 
the computed P decreases as ~ r^""*^ (see bottom panel of 
Fig. 13 . Possible origins of this discrepancy are discussed 
in Sect. EH 

Another consequence of a non-gaussian dependence of 
Xk (or Xk) with the frequency is a larger relative contri- 
bution of the Reynolds stress Pr to the mode excitation 
rate P. This is shown in Fig. 0] The GF generates a rela- 
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Fig. 3. Top: The curves correspond to computed values 
of P(y) obtained with different analytical functions for 
Xk{y): the GF (long dashed curve), the GEF (dashed 
curve), the GLF (dots-dashed curve) and the LF (solid 
curve). In all calculations, we assume A = 1. The dots rep- 
resent P(y) derived from the amplitudes and line widths 
of the f = p modes measured by ChaDlin_ct„aI. (1998). 
Bottom: same as the top panel but P is plotted in a 
log-log representation as it is usually represented in the 
literature. The vertical and horizontal scales have been 
c hosen for a easy comparis on with equivalent plots found 
in'Stein , & Nordlundl 1 20nitl . The lines with dots show two 
different power laws : one with p = —6.2 and the other 
with p = — 1. 

tive contribution of the Reynolds stress to the excitation 
{Pfi/P) which is smaller than that obtained assuming a 
non-gaussian function (e.g. for the LF, the relative contri- 
bution of to the excitation is at least ^ 2 times larger 
than with the GF). Excitation by the entropy fluctuations 
takes place predominantly at the top of the excitation re- 
gion over a thin layer (< 0.2 Mm) while that due to the 
Reynolds stress extends deeper below (~ 0.5 — 2 Mm). At 
the top of the excitation region, the discrepancy between 
the GF and Xk inferred from the 3D simulation mainly 
occurs above ~ 5 mHz and thus has a small impact 



Fig. 4. Top: Same as Fig.Ofor the relative contribution 
of the Reynolds stress, Pr to the total acoustic energy P. 

on mode excitation. This is not the case deeper in the 
excitation region where the GF under- estimates Xfe in ^ 
frequency range increasing inward. 

5. Conclusion and discussion 

5.1. A Non-Gaussian eddy time correlation 

In the present work we characterize empirically Xk and of 
X%, the frequency components of the correlation product 
of the turbulent velocity field and of its vertical component 
respectively. A frequency analysis of a solar 3D simulation 
shows that at large scales (fc ~ ko) the gaussian function 
significantly underestimates Xk and Xk frequency 
range (i^ ~ 2 — 4 mHz) where acoustic energy injected 
into the solar p modes is the largest. 

As a result, the maximum value of P is found ^2.7 
smaller than the solar seismic constraints. 

This partly explains the underestimate of the values of 
solar p mode excitation rates obtained by iHoudek et alJ 
lll999tl whose coi nputations are b ased on the theoreti- 
cal expression by iBalmforthI l)l992|) and by ISamadi et alJ 
(2003). 

In order to reproduce the main properties of Xk 
(or Xk) ' one has to consider a model which must 
decrease more slowly with i' than the GF. We then 
assume for Xk and xt three different simple analyti- 
cal forms : the Lorentzian Function (LF), the so-called 
'Gaussian Exponential Function' (GEF, which is com- 
posed by the GF plus an exponential function) and the 
so-called 'Gaussian Lorentzian Function' (GLF, which is 
composed by the GF plus a Lorentzian function). 

From the top of the excitation region (which corre- 
sponds to the top of the superadiabatic region) down to 
the middle of the excitation region (z ^ —0.5 Mm where z 
is the distance to the radius at the photosphere), the best 
agreement between Xk and the analytical approximations 
is obtained with the LF and with A = 1. Deeper within 
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the excitation region (z < —0.5 Mm), the agreement is 
better for A < 1 and A decreasing with depth. 

The frequency dependencies of xl ^^'^ Xk are found to 
be very similar. However xt is best modeled by the GEF. 
As for Xk, the agreement is better below z ~ —0.5 Mm 
with decreasing values of the parameter A than with A = 1. 

Assuming a non-gaussian function - cither the LF, the 
GEF or the GLF - results in values for Pmax, the maxi- 
mum of excitation power, which are 2 times larger than 
when assuming the GF and brings Pmax much closer to 
the maximum of P d erived from the solar seismic data of 
iGhaDhn et all lll998^ . 

We also find that taking into account the variation of 
A with depth for z below —0.5 Mm does not significantly 
change the values of P. A constant value can then be 
assumed in the calculation of the solar p mode excitation 
rates. The constant value of A on the other hand plays an 
important role and we find A = 1. 

We have investigated the sensitivity to the adopted 
representation for Xk ■ Although the LF fits best the ly- 
variation of Xk inferred from the 3D simulation, the GLF 
results in value for Pmax closer to the seismic constraints. 
However, the differences obtained with the different non- 
gaussian approximations for Xk are globally smaller than 
the actual err o r bars associated with the observations of 
iGhaphn et al.l l(l998t) . On the other hand, below the fre- 
quency range where observational constraints on P arc 
available (i.e. below v < 1.8 mHz), the differences be- 
tween P obtained with different non-gaussian functions 
are very large compared to the current error bars. Those 
differences are directly related to the diffences in the ^- 
variation of the non-gaussian forms investigated in this 
work. This suggests that accurate enough data below this 
frequency range, could provide confirmation that the LF 
is indeed the best model for Xk- 

5.2. Relative contribution of tfie entropy fluctuations 
to the excitation 

The non-gaussian character of Xk causes the excitation 
region to extend deeper (~ 500 km for modes of or- 
der n = 20) than with the GF 200 km resp.). The 
largest entropy fluctuations mainly occur at the outermost 
part of the convective zone (CZ) over a very thin region 
(~ 100 km) while excitation by the Reynolds stress con- 
tribution occurs on a more extended region. Consequently 
the non-gaussian property of Xk leads to a relatively 
larger contribution of the Reynolds stress to the excita- 
tion than in the case of a GF. As a result, the Reynolds 
stress contribution is of the same order as the contribu- 
tion arising from the advection of the turbulent fluctu- 
ations of entropy by the turbulent movements (the so- 
called entrop y source term). This is in contrast with pre- 
vious results ijSamadi et al.ll200l|) based on the GF which 
concluded that the entropy source term dominates the 
Reynolds stress by about ~ 20. It also differs with results 
bv lGoldreich et al.. (,1994) who found that the excitation 



arising directly from the entropy fluctuations dominates 

by about ~ 10. 

( >! the other hand, in IStein fc Nordlundl ()200l|) . the 
excitation by turbulent pressure (Reynolds stress) is found 
dominant (~ 4 times larger) whereas here we find that 
the contribution of the entropy source term cannot be ne- 
glected. Whether this is the signature of some deficiency 
in the present excitation model is an open question. 

5.3. Summary 

We show that the usually adopted gaussian function for 
Xk is neither consistent with the properties of Xk inferred 
from the 3D simulation nor does it reproduce the observed 
maximum of the solar p-modes excitation rates. 

Following an empirical approach we improve the model 
of the convective eddy time-correlation Xk which enters 
the current model of stochastic excitation. We then show 
that to reproduce both the z^-variation of Xk as inferred 
from the 3D simulation and the observed maximum of the 
solar p-modes excitation rates one has to consider a non- 
gaussian form which decreases at high frequency slower 
than the GF, as do the different non-gaussian functions 
investigated here. 

The use of non-gaussian functions, for instance the 
LF, reproduces reasonably well the maximum value of 
the rate at which solar p-modes are excited without any 
ajdustements of free parameters or without introducing 
a sc aling factor, in contrast with previous approaches 
re.g.lBalmforthlll992l:lc^ldreich et al.lll994l:ISamadi et alJ 
l200l|) . We then solve the problem of the underestimation 
by the previous theoretical approaches. Furthermore the 
use of such a non-gaussian form for Xk makes the contri- 
bution of the turbulent pressure to the excitation much 
larger than in previou s works making our resu lts more 
consistent with that bv lStein fc Nordlundl l)200l|) . 

Our investigation clearly emphasizes the non-gaussian 
character of the solar p modes excitation as a result of 
the non-gaussian property of the convective eddies time- 
correlations. It also shows that the dynamic properties of 
the solar turbulent convection inferred from the 3D simu- 
lation are consistent with the helioseismic data. 

We stress that only simple non-gaussian forms for Xk 
have been investigated here. More sophisticated forms are 
likely to improve the agreement with the z^-dependency 
of Xk (or Xk)- This would not affect the main conclusions 
presented in the present paper. 

5.4. Possible origin of the non-gaussian property of Xk 

We recall that Xk measures the temporal evolution of the 
correlation between two points of the turbulent medium 
separated by a distance of ~ 27r/fc. A gaussian time- 
correlation means that the fluid motions in the medium 
are random in time. Departure from a gaussian time- 
correlation at large scales (fc ~ fco) suggests that a strong 
correlation exists at that scale. 
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Downward plumes are likely to be responsible for 
the non-gaussian behaviour of Xk- Downward and up- 
ward convective motion s are indeed highly asymmetric 
llStein & Nordhmdlll998l) : downward flows are associated 
with patterns (plumes) whi ch are more coherent tha n the 
upward moving structures ijRieutord fc ZahiJll995|) . The 
upward flows are associated with less coherent and more 
random structures (granule s) characterised by a bro ad va- 
riety of sizes and lifetimes l)Rieutord fc Zahnlll995() . The 
non-gaussian behavior of Xfc can most probably be at- 
tributed to plumes. This however remains to be checked 
(work in progress). 

5.5. Possible origin of the remaining discrepancy 

Despite a clear improvement in the agreement between 
observed and theoretical excitation rates, important dis- 
crepancies between the computed P and the solar mea- 
surements still remain at high frequency v > 3.5 mHz 
(see Sect. Eland Fig.OJi. 

On the 'observational side', at high frequency, larger 
uncertainties for the damping rates r) induce larger uncer- 
tainties on the derived supply energy rates. 

On the theoretical side, part of the discrepancy might 
well be attributed to a poor description of the eigenfunc- 
tions at high frequency. Indeed, the discrepancies between 
the calculated eigenfrequencies and the observed ones are 
largest at high frequency (i^ > 3 mHz). This indicates 
that the description of the eigenfunctions are less accu- 
rate at high frequency. As the expression for calculating 
P involves the first and second derivatives of the mode 
eigenfunction, the lack of accuracy in the calculation of 
the eigenfunctions has a larger impact on P at high fre- 
quencies than at small frequencies. 

Other possible causes can perhaps be related to our 
simplified excitation model which assumes isotropic turbu- 
lence. Indeed the current theory assumes that the stochas- 
tic excitation is the same in all three directions, par- 
ticularly between the ascending and descending flows. 
However the kinetic energy and entropy fluctuations are 
larger in the downward flows than in the upward flows 
l)Stein fc Nordhmdlll99j^ . Therefore the driving arising 
from the advection of the turbulent fluctuations of entropy 
by the turbulent movements differs significantly between 
the elements moving downwards and those moving up- 
wards. As the entropy fluctuations are largest in the out- 
ermost part of the convective zone, the above mentioned 
asymmetry will predominantly affect the high frequency 
modes. 

Moreover, it is also assumed that the total kinetic 
energy, E, is isotropically injected in all 3 directions. 
Excitation of the radial p modes results from the ver- 
tical component of the velocity. However at the top of 
the convective zone, the distribution of kinetic energy in 
E[k, z) and in Ez{k, z) are very different from each other. 
These differences may affect more strongly the high fre- 



quency modes. Consequences of these departures from the 
isotropic assumption need to be further investigated. 

5.6. Perspectives 

The non-gaussian property of Xk and its consequences for 
the stochastic excitation has been investigated so far only 
for the Sun. However such a non-gaussian feature of the 
turbulence will most likely also be of importance for solar- 
like oscillating stars more massive than the Sun, provided 
our analysis is also valid for these stars. This can sub- 
stantially change the excitation spectrum P for such stars 
compared to that which is currently predicted. 

Therefore investigations of p mode excitation in hotter 
and more massive stars must be undertaken, which should 
proceed in two steps : first, the validity of the present re- 
sults obtained in the solar case must be investigated for 
other stars with, for instance, the help of dedicated 3D 
simulations. The conclusions which will drawn from this 
first step must be used in a second step to study the fre- 
quency dependence and the magnitude of P for differ- 
ent solar-like oscill ating stars (see preliminary results in 
ISamadi et ai]l2n02l) . 

Future space missi ons such as COROT 

jiBaglin fc The Corot Tea^ Il998ll. MOST ijMatthewd 
1993) and Eddington l|Favata et al 1 l2nn(t will provide 
high-quality data on seismic observations. COROT will 
be the first mission that will provide high precision mode 
amplitudes and linewidths in other stars. This high- 
quality data will allow us to derive the excitation rate P. 
and will provide improved observational constraints on 
the theory of stochastic excitation which is, at present, 
poorly constrained by observation. 
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